dbc90b24-73d1-4f92-8141-44c44321c9da 35 4 5 decimal

&[DATE] &[TIME] - &[FILENAME]

Part 1: otherwise dressed from 'Carlos'=====================

N 6 :

i 1 N range f i el x i 1 -^: for

m n dot x f x g * x 1 - 1 int ≡

f transpose 1 x x 2^x 3^x 4^x 5^16 mat

o 1 el f 1 el :

<= pre-initialise

kernel algo: MUST BE pre-initialised o i 2 f rows range d 0 : j 1 i 1 - range a f i el o j el * x 1 - 1 int o j el o j el * x 1 - 1 int / 15 round : d d a o j el * - : 21 line for o i el f i el d + (: 31 line for :

Explain:

1. "Syntax incorrect" at the user level but the system of 'o' is evaluated at the kernel level, awaiting the indexed call(s) of oi. 2. The ratio of the two dot product needs be vectorised, Smath does not do it from normal code. The 'round(_,15)' has same vectorising effect, thus enabling the evaluation of the 'a(s)' coefficients, which coefficients are pure numeric. 3. The pest remains about extracting meaningful components for export and more usefull sub-use. However OK for QuickPlot. 4. Lagrange "roots style" export is compatible for Engineering works, mostly fitting sampled data. Think in term of Fourier spectrum from clean model.

Collect/Normalise function style

x o o :

x o 1 el 1

x o 2 el x

x o 3 el 500000000000000 x 2^* 166666666666667 - 500000000000000 /

x o 4 el x 150000001 - 250000000 x 2^* + (* 250000000 /

x o 5 el 5568434188608080100000000666667 - 500000000000000 x 4^* + (* 487229337000119166666666666667 - 500000000000000 x 2^* + (* - (* 1421085471520200000000000000000 /

Poly x o 1 el x o 2 el x o 3 el x o 4 el x o 5 el x o 6 el 61 sys :

Poly

The ROBUST root(s) Solve scanner f x0 x1 TOL Solve ρ x f x x0 x1 solve : i 1 ρ rows range a ρ i el 0.1 ρ i el * - : b ρ i el 0.1 ρ i el * + : xn b a + (2 / : ε TOL : n 0 : xn f abs ε > xn f a f * 0 ≤ b xn : a xn : if xn b a + (2 / : n n 1 + : 31 line while 61 line X i el xn : 21 line for X 31 line :

Export Lagrange roots style

x f x o 6 el :

TOL 1015 -^:

r f 1 - 1 TOL Solve : 0.90618 - 0.538469 - 0.538469 0.90618 41 mat

x Lagrange x r 1 el - (x r 2 el - (* x r 3 el - (* x r 4 el - (* x r 5 el - (* :

x Lagrange

Part 2: dressed tutor & export style 'Jean'===================

FUNCTION

m n dot m n * x 1 - 1 int m m * x 1 - 1 int / 15 round :

The method exposed here is piece wise, more effort demanding but convenient for Engineering works ... visible.

f transpose 1 x x 2^x 3^x 4^x 5^16 mat

o1 f 1 el :

o1 1

a2 o1 f 2 el dot :

a2 0

o2 f 2 el a2 o1 * - :

o2 x :

a3 o1 f 3 el dot : 0.333333

b3 o2 f 3 el dot : 0

a3 0.3333

b3 0

o3 f 3 el a3 o1 * - b3 o2 * - :

o3 x 2^13 / - :

a4 o1 f 4 el dot : 0

b4 o2 f 4 el dot : 0.6

c4 o3 f 4 el dot : 0

a4 0

b4 0.6

c4 0

o4 f 4 el a4 o1 * - b4 o2 * - c4 o3 * - :

o4 x 3^35 / x * - :

a5 o1 f 5 el dot : 0.2

b5 o2 f 5 el dot : 0

c5 o3 f 5 el dot : 0.857143

d5 o4 f 5 el dot : 0

a5 0.2

b5 0

c5 0.857143

d5 0

o5 f 5 el a5 o1 * - b5 o2 * - c5 o3 * - d5 o4 * - :

o5 x 4^335 / + 67 / x 2^* - :

combine the terms 67 / 13 / - (* - 15 / - 335 /

5 first orthogonal functions

o1 o2 o3 o4 o5 5 1 sys

NOT convenient symbolic o1 0203 o4 0551 sys 10203 x 3 - 5 x 2^* + (* 5 / 0551 sys

o1 o2 03040551 sys 1 x x 2^13 / - x 3^35 / x * - x 4^335 / + 67 / x 2^* - 51 mat ≡

x o1 x f1 :

x o2 x f2 o1 f2 dot o1 o1 dot / x o1 * - :

x o3 x f3 o1 f3 dot o1 o1 dot / x o1 * - o2 f3 dot o2 o2 dot / x o2 * - :

x o4 x f4 o1 f4 dot o1 o1 dot / x o1 * - o2 f4 dot o2 o2 dot / x o2 * - o3 f4 dot o3 o3 dot / x o3 * - :

x o5 x f5 o1 f5 dot o1 o1 dot / x o1 * - o2 f5 dot o2 o2 dot / x o2 * - o3 f5 dot o3 o3 dot / x o3 * - o4 f5 dot o4 o4 dot / x o4 * - :